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Abstract.  We  introduce  a  new  class  of  exponentially  accurate  filters  for  processing  piecewise 
smooth  spectral  data.  Our  study  is  based  on  careful  error  decompositions,  focusing  on  a  rather 
precise  balance  between  physical  space  localization  and  the  usual  moments  condition.  Exponential 
convergence  is  recovered  by  optimizing  the  order  of  the  filter  as  an  adaptive  function  of  both  the 
projection  order,  and  the  distance  to  the  nearest  discontinuity.  Combined  with  the  automated  edge 
detection  methods,  e.g.,  [GeTa02],  adaptive  filters  provide  a  robust,  computationally  efficient,  black 
box  procedure  for  the  exponentially  accurate  reconstruction  of  a  piecewise  smooth  function  from 
its  spectral  information. 
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1.  Introduction 

The  Fourier  projection  of  a  27r-periodic  function,  Sj\rf(-),  enjoys  the  well  known  spectral  convergence 
rate,  that  is,  the  convergence  rate  is  as  rapid  as  the  global  smoothness  of  /(•)  permits.  Specifically,  if 
/(•)  has  s  bounded  derivatives  then  \S^f(x)  —  f(x) \  <  Const ||/||cs  •  IV1'5.  This  interplay  between 
global  smoothness  and  spectral  convergence  is  reflected  in  the  dual  Fourier  space  through  the 
rapidly  decaying  Fourier  coefficients  \f(k)\  <  2Tr\\f\\cs\k\~s .  On  the  other  hand,  spectral  projections 
of  piecewise  smooth  functions  suffer  from  the  well  known  Gibbs’  phenomena,  where  the  uniform 
convergence  of  S]yf(x)  is  lost  in  the  neighborhood  of  discontinuities  and  the  global  convergence  rate 
of  Sj\rf(x)  deteriorates  to  first  order.  These  related  phenomena  are  manifestations  of  unacceptably 
slowly  decaying  Fourier  coefficients. 

Two  interchangeable  processes  are  available  for  recovering  the  rapid  convergence  in  the  piecewise 
smooth  case.  These  are  mollification,  carried  out  in  the  physical  space  and  filtering,  carried  out  in 
the  Fourier  space.  Filtering  accelerates  convergence  when  premultiplying  the  Fourier  coefficients 
f{k)  by  a  rapidly  decreasing  function  cr (• ) ,  resulting  in  modified  coefficients,  f(k)a(\k\/N),  with  a 
greatly  accelerated  decay  rate  as  \k\  |  A.  This  accelerated  decay  in  the  dual  space  corresponds  to  a 
smoothly  localized  mollification  in  the  physical  space.  In  [TT02]  we  showed  how  to  parameterize  an 
optimal  mollifier  in  order  to  gain  the  exponential  convergence  for  piecewise  analytic  /’ s.  The  key 
ingredient  in  our  approach  was  adaptivity,  where  the  optimal  mollifier  is  adapted  to  the  maximal 
region  of  local  smoothness.  Here  we  continue  the  same  line  of  thought  by  introducing  adaptive 
filters,  which  allow  the  same  optimal  recovery  of  piecewise  smooth  functions  from  their  Fourier 
coefficients.  In  particular,  piecewise  analytic  functions  are  recovered  with  exponential  accuracy.  A 
brief  overview  follows. 

We  consider  a  family  of  general  filters  u(-)  which  are  characterized  by  two  main  properties.  First, 
we  seek  the  rapid  decay  of  ay  :=  a(\k\/N)  which  is  tied  to  a  regular,  compactly  supported  multiplier 
a  £  Cq[—  1, 1],  Being  compactly  supported,  such  filters  are  restricted  to  A-Fourier  expansions, 
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The  operation  of  such  filters  in  Fourier  space  corresponds  to  mollification  in  physical  space,  ex¬ 
pressed  in  terms  of  the  associated  modifier,  <hcr(y)  :=  1/2-k  cr(\k\/N)elky , 

(1.2)  sft/(x)  =  f* ^(x)  =  Y  f  - y)dy  *a(y)  ■=  j-  E 

71  71  |fc|<JV 

Second,  such  filters  are  required  to  satisfy  the  usual  moments  condition,  e.g.  [MM078,  Va91] 

(1.3)  f  yn<$>a{y)  =  6n0,  n  =  0, 1, . . .  ,p  -  1  <  q. 

J  —  TV 

The  first  requirement  of  CQ-smoothness  is  responsible  for  localization  —  the  essential  part  of  the 
associated  modifier,  is  supported  near  the  origin,  consult  (2.5)  below.  The  second  property 
drives  the  accuracy  of  the  filter,  by  annihilating  an  increasing  number  of  its  moments. 

The  rich  subject  of  filters  include  the  classical  filters  of  finite  order  accuracy  where  finite  p  <  q 
dictate  a  fixed  convergence  of  polynomial  order,  0(N~P),  consult  [Va91].  By  letting  q  f  oo  one 
obtains  a  Cq°[— 1,  l]-filter,  that  is,  an  infinitely  differentiable  compactly  supported  filter  a,  which 
respects  (1.3)  for  increasing  orders  p.  Majda  et.  al.,  [MM078]  employed  such  filters  to  postprocess 
piecewise  solutions  with  propagating  singularities  and  achieve  spectral  convergence  in  the  sense 
of  having  a  convergence  rate  faster  than  any  fixed  order.  Vandeven,  [Va91],  constructs  spectrally 
accurate  filters  by  relating  the  order  of  the  filter,  ( q,p ),  to  the  increasing  order  of  the  projection, 
q  =  q(N),p  =  p(N).  An  alternative  approach  for  spectral  accuracy  employs  highly  oscillatory 
modifiers  which  are  activated  in  physical  space.  In  [GoTa85],  Gottlieb  and  Tadmor  constructed 
such  (properly  dilated)  modifiers  of  the  form  <d(y)  =  p(y)Dp(y),  where  Dp(-)  stands  for  the  usual 
Dirichlet  kernel  of  degree  p  ~  y/N  and  p  is  a  standard  Cq°[—  1, 1]  cut-off  function  normalized  such 
that  p(0)  =  1. 

The  different  filters  and  modifiers  advocated  in  these  works  enable  to  reconstruct  the  underlying 
piecewise  smooth  data  from  its  given  spectral  content.  Spectral  accuracy  is  achieved  in  smooth 
regions  as  long  as  they  are  bounded  away  from  discontinuities,  but  the  error  deteriorates  in  the 
neighborhood  of  such  discontinuities  due  to  spurious  oscillations.  The  latter  difficulty  was  addressed 
by  Gottlieb,  Shu  and  collaborators,  by  invoking  Gegenbauer  expansions  which  are  driven  by  a 
judicious  choice  of  a  localizer  (1  —  y2)x  which  are  appended  to  the  Dirichlet  kernel  Dp(y ),  consult 
[GoSh98]  and  the  references  therein.  Their  approach  allows  for  high  resolution  uniformly  up  to  the 
discontinuities,  but  its  precise  (p.  A)-parameterization  as  a  function  of  N  has  a  rather  sensitive  /- 
dependence  which  impact  the  overall  robustness  of  the  Gegenbauer  reconstruction,  e.g.,  [Bo05].  In 
[TT02]  we  have  introduced  an  alternative  approach  where  the  accuracy  is  adapted  according  to  the 
maximal  region  of  local  smoothness.  Specifically,  we  have  shown  how  the  Gottlieb-Tadmor  mollifiers 
are  optimized  when  their  order  is  chosen  adaptively  as  p  ~  Nd(x).  Here  d{x)  is  distance  between  the 
location  x  to  its  nearest  discontinuity,  d(x)  =  distance(x,  sinsuppf  (•));  the  distance  function  d(x) 
could  be  recovered  from  the  Fourier  coefficients  by  edge  detection,  e.g.,  [GeTa00],[GeTa02],  The 
resulting  adaptive  mollifiers  lead  to  exponentially  accurate,  numerically  robust  mollifiers  of  order 
exp(—a(KNd(x))1^°l)  with  a  >  1  dictated  by  the  detailed  Cq° -regularity  of  p;  specifically,  a  >  1 
reflects  the  Gevrey  regularity  of  p  (for  Gevrey  regularity  and  the  similar  class  of  ultramodulation 
spaces  we  refer  to  e.g.,  [Jo]  and  [PT02],  respectively).  The  key  ingredient  in  our  adaptive  approach 
is  giving  up  the  exact  moments  condition;  instead,  (1.3)  is  satisfied  modulo  these  exponentially- 
negligible  errors,  replacing  the  exact  (1.3)  with  the  requirement 

(1.4)  a(n)(0)  =  dn0,  n  =  0,1,...  ,p  —  1  <  q. 

The  precise  relation  between  (1.4)  and  (1.3)  is  quantified  in  theorem  2.2  below.  We  note  that  it 
is  rather  simple  to  construct  admissible  filters  satisfying  the  last  requirement  for  an  arbitrary  p;  a 
prototype  example  is  given  by  the  Cg°[—  1,  l]-filters 
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(1.5)  *P(Q  =  (  (f£l)  ’  lfl  <  1 

lo  iei>i- 

The  purpose  of  this  paper  is  to  construct  a  new  class  of  exponentially  accurate  adaptive  filters. 
As  before,  the  key  issue  is  the  parameterization  of  their  order,  p.  Here  we  develop  the  rigorous 
study  for  the  optimal  parameterization  for  such  filters.  We  advocate  adaptive  filters  in  the  sense 
that  their  order,  p  =  p(N,d(x)),  depends  on  both  —  the  order  of  the  projection,  N,  and  the 
distance  function  d{x).  Summarized  in  theorem  2.1  below,  our  main  result  states  that  the  optimal 
adaptive  filter  is  determined  to  be  of  order  p(N)  ~  (Nd(x))1^a  with  a  >  1  reflecting  the  Gevrey 
regularity  of  a.  While  achieving  exponential  accuracy  away  from  discontinuities,  the  new  filters 
are  adapted  so  as  to  prevent  spurious  oscillations  throughout  the  computational  domain,  including 
discontinuous  neighborhoods.  We  mention  here  the  adaptive  filters  introduced  by  Boyd  in  [Bo96]. 
Boyd’s  procedure  was  based  on  the  acceleration  summability  by  the  so  called  Euler  lag  averaging; 
the  acceleration  was  limited,  however,  since  the  resulting  piecewise  constant  filters  of  order  p  ~ 
Nd(x)  were  consistently  larger  than  the  optimal  order  and  they  exhibit  slower  convergence  than 
the  non-adaptive  order  of  [MM078]. 

Our  current  discussion  on  adaptive  filters  follows  a  similar  approach  for  the  adaptive  Gottlieb- 
Tadrnor  mollifiers  constructed  in  [TT02],  p(y)Dp(y),  where  the  precise  Gevrey  regularity  of  the  p 
allows  us  to  obtain  tight  error  bounds  which  in  turn  reveal  the  optimal  adaptive  parameteriza¬ 
tion,  p  =  p(N,d(x)).  New  tight  error  bounds  are  outlined  in  §2  and  are  confirmed  by  numerical 
simulations  in  S3. 


2.  Adaptive  order  filters 

In  this  section  we  show  how  the  regularity  and  moments  properties  of  the  filter  a  are  translated 
into  precise  statements  of  localization  and  accuracy  of  the  associated  modifier  &a(x).  We  begin  by 
decomposing  the  filtering  error  /(•)  —  =  /  —  /  *  <h<T  into  the  two  terms 


(2.1) 


f(x)  -  f  * 


[  $a(y) L/0)  -  f(x  -  y)}  •  [i  -  x(y)]dy  + 

J  —TV 

+  [  $CT(y)[/0r)  -  f(x  -  y)\x(y)dy  =■■  2i  +t2. 

J  —TV 


Here,  %(•)  =  Xx(')  is  a  auxiliary  cut-off  function  adapted  to  the  smoothness  region  of  /.  To  this 
end  we  let  d(x)  denote  the  distance  between  x  and  its  nearest  discontinuity  so  that  the  y-function 
f(x)  —  f{x  —  y)  remains  smooth  for  the  largest  symmetric  interval,  \y\  <  d[x).  We  then  set 
x(y)  =  Xx(y)  :=  p(y/d(x))  where  p  is  a  standard  Cq°  cut-off  function, 

X(y)  ^  :=  p(jjfi)  p(y)  =  {  =  ^  j*j  ;  f 

We  observe  that  the  dilated  cut-off  function  Xx{y)  =  p(y/d(x))  enforces  the  support  of  the  first 
integrand  on  the  right  of  (2.1)  to  be  bounded  d(x) /2-away  from  x  while  the  second  term  is  supported 
in  the  d(x)-neighborhood  of  x.  To  simplify  matters,  we  assume  that  p  is  adapted  to  the  same  Cq° 
regularity  of  a. 
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We  turn  to  estimate  the  first  error  term  on  the  right  of  (2.1)  which  measures  the  essential  localization 
of  the  mollifier.  To  this  end,  we  use  the  following  aliasing  formula,  expressing  our  IV-degree  modifier, 
$CT(y)  in  terms  of  the  equally  sampled  inverse  Fourier  transform,  pa{y)  :=  f  a(£)eiy^d£( 1 ), 

N  °° 

(2.2)  <f>a(y)  =  —  Y  pa({N{y  +  2nn)),  4>CT(y) 

n=— oo 

The  usual  Fourier  decay  rate  estimates  then  yields, 

t.  j  OO  CAJ 

|3>f7(?/)|  <  ^  Y  l<^'T(iV(2/  +  2™))l  <  N^WaWcp  Y  I y  +  2nn\~p 

n=— oo  n=— oo 

(2.3)  <  Const- N\\cr\\cp(N\y\)~p,  Vp,  j/f e[-7T,7t]. 

We  observe  that  the  first  integrand  on  the  right  of  (2.1)  is  supported  across  the  possible  dis¬ 
continuities  of  f(x)  —  f(x  —  •).  Lack  of  smoothness  excludes  the  possibility  of  high-oscillatory 
cancelations.  Instead  we  now  seek  a  tight  upper  bound  on  the  decay  of  the  associated  mollifier, 
for  \y\  >  d(x)/2.  To  this  we  need  to  quantify  the  C°° -regularity  of  our  filter  a.  We  focus  on  <r’s 
which  have  Gevrey  regularity  of  order  a,  denoted  Ga  below;  in  our  case,  a  =  ap  in  (1.5)  belong  to 
G2,  namely,  there  exist  constants,  M  =  Ma  and  y  =  r]a  >  0  (independent  of  p)  such  that 

(2.4)  hp\\cp  <  Ma(p\)ay~p,  a  =  2,  <tp(£)  =  (  GXP  (^2  -  1 

lo  |CI>i. 

Details  are  outlined  to  lemma  2.1  below.  Incorporating  the  above  growth  rate  into  the  localization 
estimate  (2.3)  yields 


l^(l/)|  <  Const. Ma(pl)2  ( — -j-r — )P, 

^v<TN\y\' 

which  is  minimized  at  p  =  pmin  :=  (77  ■  Nd{x))1^2 .  This  shows  that  with  this  choice  of  adaptive  p, 
the  mollifier  associated  with  our  u-filter,  is  essentially  localized  in  the  neighborhood  of  x,  as  it 
admits  an  exponential  decay 

(2.5)  |4>CTp(y)|  <  Const a  ■  Ne~ ^N^1/2 . 

Here  and  below  ij  is  a  positive  constant  which  may  differ  among  the  different  estimates.  In  particu¬ 
lar,  since  [1  —  Xx(y)]  and  hence  the  first  integrand  on  the  right  of  (2.1)  are  supported  at  |y|  >  d(x)/ 2, 
the  exponential  bound  follows 

(2.6)  |Ii  |  <  Const aj  ■  Ne~(p°Nd{x))1/2 

We  turn  to  the  second  error  term,  X2  =  f  <S>a(y)[f(x)  —  f(x  —  y)]Xx{y)dy.  Ti’aditionally,  such  a  term 
is  upper  bounded  by  {d(x))p\\f\\CP[x-d(x),x+d(x))/p-  through  Taylor  expanding  f(x  -  y)  about  y  =  0 
and  by  invoking  the  moments  condition  (1.3).  This  bound  is  useful  for  a  vanishing  neighborhood, 

1The  result  follows,  for  example,  by  sampling  the  Fourier  transform  <r(£)  =  1/27T  J  {y)e~ly^dy, 

°  (  jy)  =  £/  (Ny)e-%Ny^  dy 

=  5“  f  ^ {N (y  +  ^n))e~lWy dy  = I  (  Y  ^ (N {y  +  2nn))j  e~Ak{y dy , 

n=  —  00  ~7r  ~‘K  \n=  —  00  / 

and  comparing  with  the  discrete  inverse  Fourier  transform,  a(\k\/N)  =  f  4>CT (y)e~^k^y dy 
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d(x)  <C  1,  while  suffering  by  increasing  the  contribution  of  the  first  term  on  the  right  of  (2.1),  as 
reflected  through  its  upper-bound  (2.6).  We  therfore  let  d(x)  be  as  large  as  possible  so  we  cannot 
argue  by  localization.  Instead,  this  portion  of  the  error  decreases  due  to  cancelation  of  oscillations 
by  increasing  the  order  p  of  <&av .  To  this  end  we  write 

(2.7)  12  =  I  ^ap(y)g(y)dy  =  ^2  a  f^P\g(k),  g{y)  =  gx{y)  :=  [f(x)  -  f(x  -  y)\Xx{y), 

J~'K  |fc|<JV  '  ' 

and  we  turn  to  estimate  the  Fourier  coefficients  on  the  right.  By  our  assumption,  f{x)  —  f(x  —  y) 
remains  analytic  for  |y|  <  d(x)  and  hence  gx(y)  =  [f(x)  —  f(x  —  y)]x(y)  is  C°° .  We  quantify  the 
C°° -regularity  of  Xx(')  in  terms  of  the  same  Gevrey  regularity  of  order  a  =  2  that  a  has,  so  that 
HXxllo  <  M(p\)2(r)pd(x))~p.  Thanks  to  the  analyticity  of  f{x)  —  f(x  —  y),  it  follows  that  if  p(-) 
and  hence  Xx(')  belong  to  Gevrey  class  Ga,  so  does  gx(y)  =  [/(x)  —  f(x  —  y)]Xx{y)i  and  hence 

\\9x(y)\\c*  <  M ^[x)g)P'  12/1  <  d(x)' 

The  constants  M  =  Mp  r)  and  q  =  rjpj  capture  the  detailed  Gevrey  and  analyticity  properties  of 
p(y)  and  f(x  —  y)  for  \y\  <  d(x)\  the  order  p  is  arbitrary.  The  Fourier  coefficients  g(k)  in  (2.7)  do 
not  exceed 


( 7)1)2 

\9(k)\  <  Const. \\gx(y)\\Cp\k\-p  <  Const.M  ,  g  =  gpJ. 

For  a(\k\/N),  we  distinguish  between  the  low  modes,  \k\  <  N/2  and  the  high  modes  N/2  <  \k\  <  N, 
setting 


12\  ■■= 


^22  '■  = 


E  h(v 

| fc  |  <  JV/2  V 


E 

N/2<\k\<N 


—  1 


1  -  )  -  ! 


g(k) 
g{k) 


Since  g(y)  =  [/(•)  —  /(•  —  y)]x(y)  vanishes  at  y  =  0  we  have  ~^2g{k)  =  f/(0)  =  0  and  hence 
12  =  ^21  +2^22  T  2^23  where  I2'.i  '■=  —  ^2\k\>N^(k)-  For  the  first  term  we  use  a  Taylor  expansion 
around  the  origin:  the  accuracy  assumption  (1.4)  yields 


CTrn 


< 1  ii  ii  my 


Restricted  to  the  [—1/2, 1/2]  interval,  a  retains  an  analytic  bound  ||<r||crp  <  Const.p\r]ap,  and  hence 
2T2i  does  not  exceed 


|2”21 


(2.8)  < 


E  [■ 

\k\<N/2 


-  1 


g{k) 


< 


O' 


pi 


cni-li  d 


E 

| fe|  <  JV/2 


Const.N{p\y 


(r]Nd(x))P  ’ 


V  =  VaVpj- 


(my  (P\)2 

\N  )  (' VP,f\k\d{x))P 


For  the  high  modes,  'gx(k)  are  sufficiently  small  so  that  the  simple  bound  of  |<tp(| fc | /2V) |  <  1  will  do 
for  X22, 
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(2.9) 

\M  ■■= 


E 

N/2<\k\<N 


(Jr) 


g(k) 


< 


E 


(p')‘ 


N/2<\k\<N 


(pa\k\d{x))P 


<  Const.N(p\y 


1 


(gCTNd(x)/2)p' 


Similarly,  X23  does  not  exceed 


(2.10) 


^23 1  := 


E  9(k) 

|fc|  >JV 


<  (p!)2  E 


1 


\k\>N 


(' r)a\k\d(x))P 


<  Const. N(p\y 


1 


(rjaNd(x))P ' 


We  combine  the  last  three  bounds  to  conclude 


|J2|  <  Const. N{p\)2  j  ,  g  =  mm(riar)pj,r)a/2,r)pj) 

which  is  minimized  at  the  same  value  as  before,  p  =  pmin  '■=  (rl '  Nd{x))l( 2 ,  so  that 

(2.11)  |X2|  <  Const  •  Ne~(llNd{x))1/2 

Finally,  we  recall  the  assumed  regularity  of  p  is  in  fact  dictated  by  that  of  cr  and  hence  the  various 
bounds,  p  =  Po.f  ■  We  summarize  by  stating 

Theorem  2.1.  Given  the  Fourier  projection  Sjyf  of  a  piecewise  analytic  function  /(•),  we  consider 
a  [-1-1]  filter  ct(£), 


Assume  that  a  has  G a-regularity  and  that  it  is  accurate  of  order  p  in  the  sense  of  satisfying  the 
moments  condition 


(2.12) 


crdd(o)  =  Sno,  n  =  0, 1, . . .  ,p  -  1. 


We  set  the  adaptive  order  p(x)  :=  (p-Nd(x))1(a  depending  the  distance  function  d{x)  =  dist(x,  sinsupp f  (•)) . 
The  resulting  adaptive  filter,  Sfjf,  recovers  the  pointvalues  f(x)  with  the  following  exponential  ac¬ 
curacy 


(2.13)  |/(x)  -  SaNf(x) |  <  Const  •  Ne-^m^))1/a . 

The  constant  p  =  paj  is  dictated  by  the  specific  Gevrey  and  piecewise- analyticity  properties  of  a 
and  f. 


We  close  this  section  with  the  promised  statements  on  the  exponential  error  bound  (2.4). 

{exp  { — — —  |  p  <^2. 

y  —  1 J  ’  s  Then  there  exist 

0  iei>i- 

constants  g  such  that 


\Wp\\cp  <  Const. (p\)2 g  p, 
CJpllcp([--  -])  —  Const.p\p~p. 


(2.14) 

(2.15) 
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To  verify  (2.14)  we  first  note  that  is  a  collection  of  polynomial  terms  which  premultiply  the 
exponential  in  the  variable  £p/(£2  —  1).  Each  derivative  of  ap  doubles  the  number  of  such  terms; 

thus,  by  successive  application  of  Leibniz’s  rule,  op  consists  of  2s  a-terms,  each  of  which  is  of  the 
form 


(2T6)  Ca  JJ  (  *  ) (  j)exp(- ^-y) ,  \a\  =  J2  aj  =  s- 

Here,  the  Ca’s  are  constant  integers  with  \Ca\  <  for  some  fixed  r/i  >  0.  We  consider  the 
prototype  term  TP)S  :=  corresponding  to  a  =  (0,  0, . . .  ,  s), 


Tp,s 

(2.17) 


(  \w 

V£2  - l) 


exp 


p\  rv-s+k  M  (  £P  \ 

(p-a+fc/  (e  -  D^exp\e~i) 


+ 


lower  order  terms  . 


Here,  by  ‘lower  order  terms’  we  refer  to  the  singular  behavior  of  (£2  —  1  )“■?',  j  <  k  near  £  =  ±1,  which 
weaker  then  the  leading  term  (£2  —  l)-fc+1.  To  control  the  amplitude  of  Tp^s  we  let  a(£)  :=  (£2  —  1) 
and  note  that  the  expression  |a(£)|~fcexp(aa(£)  +  /?/a(£))  is  maximized  at  £  =  £max  such  that 
a{£max)  ~  —/3/A;,  yielding 


|TPiS|  <  Const.  YY 


(p  —  s  +  k)l 


k\kke  k  <  Const.pl  ^ 


k\  <  Const.2sp\s\. 


The  other  2s  terms  in  (2.16)  admit  similar  bounds  and  the  resulting  Tp.p  bound  yields  (2.14)  with 
rj  =  4r/i.  To  prove  (2.15),  we  restrict  attention  to  a  subinterval  which  is  bounded  away  from  ±1, 
so  that  the  ^-dependent  terms  in  (2.17)  remain  uniformly  bounded,  (£2  —  l)_:,exp(£p/(£2  —  1))  < 
P2k i  j  <  k  +  1,  and  we  are  left  the  desired  upper  bound 


Tpp |  <  Const.  '^2  f  ^  'j  k  <  Const. 2pp\{ji2 )  p, 

k= o  ^  ' 

and  (2.15)  follows  with  7j  =  4//i  r/2-  ■ 

The  intricate  part  in  the  construction  of  such  highly-accurate  filters  or  mollifiers  is  the  further 
requirement  for  their  localization  (in  physical  space)  or  smoothness  (in  Fourier  space).  One  cannot 
increase  the  order  p  arbitrarily  without  steepening  T0-,  or  equivalently,  without  losing  smoothness 
of  <7.  The  solution  taken  here  was  to  satisfy  the  moments  condition  approximately ,  modulo  ex¬ 
ponentially  negligible  errors  while  retaining  the  desired  smoothness  properties.  We  note  that  our 
optimal  adaptive  filter  is  essentially  localized  in  the  physical  space  in  the  sense  that  the  associated 
modifier  is  exponentially  small  for  |y|  3>  1  /N,  (2.5),  see  figure  2.1.  In  contrast,  the  adaptive 
mollifiers  constructed  in  [TT02],  p(y /d(x))Dp(y /d(x))/d(x),  were  compactly  supported  in  physical 
space  (adapted  to  the  smoothness  neighborhood  of  x)  and  only  essentially  localized  in  the  dual 
Fourier  space.  The  precise  result  is  quantified  in  the  following. 

Theorem  2.2.  Consider  the  even  filter  o  with  Gevrey  regularity  Ga  satisfying  the  p-order  accuracy 
condition  (1-4),  with  p  ~  Nl/a.  Then  the  associated  mollifier,  T0-  satisfies  the  moments  condition 
(1.3)  modulo  an  exponentially  negligible  error, 
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For  proof,  we  appeal  to  (2.2) 


yn$a(y)dy=  ^  f  ynva(Ny)dy  +  ^-^2  f  ynpa(N(y  +  2im))dy  =:  Zx  +  X2. 


I  a  ~  \o/—o  /  o  r  \  a/— a  '  /  ,  f 

ly=-n  27T  27 T  ^  Jj/=_7r 


For  the  first  term  on  the  right  we  have 


[  ynVa{Ny)dy  -  f  ynpa(Ny)dy  =:  lu+ 112. 

J  V=  —  OO  ^^"«/|vl>7T 


—  J-  <■)  I  r  aj-ty  c\  I  &  t  ~  J-x  1  — 

J y=— oo  J |?/|>7r 

We  now  have Tn  =  (— iN)no(^n'>  (0)  =  5no  by  (1.4),  where  the  usual  decay  rate  |<^<T(y)|  <  Const. ||cr||cs 
M_s  yields 

|Xi2 1  <  Const.— - 1 1 a ||cs  [  yn~sdy  <  Const. (pN)1~s(s\)a ,  n<s  —  2. 

2?r  J  TV 

The  remainder  amounts  to  a  similarly  exponentially  small  term 


jvi-s  f*  1 

<  Const. -^-\\o\\Cs  j  \y\n-^—j^dy  <Const.(riN)1~s(s\)a,  n<s 

which  is  minimized  at  s  ~  (r]N)1^a  and  the  lemma  follows.  ■ 

We  note  in  passing  that  the  last  theorem  could  be  used  as  a  starting  point  for  an  alternative  proof 
of  the  main  result  stated  in  theorem  2.1. 


Figure  2.1:  The  modifier  (top)  and  its  semi-log  plot  (bottom)  with  the  modifier  defined  from  the 
filter  (3.1)  used  in  the  numerical  experiments,  with  N  =  128  and  filter  orders  p  =  4,8,  and  12  in 
(a-b),  (c-d),  and  (e-f)  respectively. 
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3.  Numerical  Experiments 
For  the  following  examples  we  utilize  the  filter 


(3-1) 


°P(0 


exp  (fr§-)  l£l  <  1 

k  o  ici  >  i 


which  has  Gevrey  regularity  of  order  a  =  2.  Its  advocated  order  is  then  optimized  at  the  adaptive 
order,  p  =  p(x)  =  nNd(x).  For  a  given  filter,  the  free  constant  cp  should  be  selected  to  enhance 
the  immediate  localization  of  <hCT(-)  by  minimizing  ||<t||c«i.  The  value  of  such  an  optimal  cp  does  not 
permit  a  closed  form  expression;  an  approximate  condition  used  in  the  numerical  examples  below 
is  cr^2l  (1/2)  =  0,  resulting  in 


_op3  9p2  +  3p  +  14 
Cp  f  4  '  9p2  +  12p  +  4 ' 


To  allow  direct  comparison  between  our  adaptive  filters  and  the  adaptive  mollifiers  advocated  in 
[TT02],  we  concern  ourselves  with  the  two  prototypes  of  piecewise  analytic  functions,  f\ (x)  and 
I2 (t)  given  below. 


(3.2) 


(3.3) 


f  sin(x/2)  xG[0,7t) 

\  —  sin(x/2)  xG[7r,27r) 


h(x) 


(2e2x  —  1  —  e7r)/(e7r  —  1)  ®e[0,7r/2) 
—  sin(2.x/3  —  vr/3)  xG[7r/2,27r) 


The  first  function,  /i(-)>  possesses  a  mild  regularity  constant  and  a  single  discontinuity  at  x  =  it; 
consequently  d(x)  =  \x  —  tt\  for  x  €  [0,  27t] .  The  second  function,  /2O),  w&s  constructed  as  a  more 
challenging  test  problem  with  a  large  gradient  to  the  left  of  the  discontinuity  at  x  =  tt/2.  Moreover, 
lacking  periodicity  /2(-)  feels  three  discontinuities  per  period; 


d(x)  =  min(\x\,  \x  —  7r/2|,  \x  —  2tt\)  x  €  [0,  2ir\. 

For  both  functions  the  exact  Fourier  coefficients,  {f(k)}k<N,  are  given  and  then  filtered  to  re¬ 
cover  the  intermediate  pointvalues  jj{y  —  \)  for  u  =  1,2,...  2N .  Graphs  (a)-(d),  use  fixed  order 
filters,  verifying  the  well  known  fact  that  higher  order  filters  gives  superior  convergence  away 
from  discontinuities  and  lower  order  filters  near  discontinuities.  Graphs  (e)-(f)  illustrate  the  su¬ 
perior  convergence  for  the  adaptive  filter  described  in  theorem  2.1,  computed  with  adaptive  order 
p  =  p(x)  =  max(2,  l(AId(x))1//2).  We  note  in  passing  that  the  same  filter  order  is  used  for  both  /i(-) 
and  _/*2  ( • )  5  ignoring  the  different  analyticity  properties  of  f\  and  fi  (reffected  by  different  analyticity 
constants  r]f),  and  achieving  exponential  accuracy  in  both  instances.  Results  of  the  adaptive  filter 
are  contrasted  with  the  spectrally  accurate  filter  of  Vandeven,  [Va91],  where  the  variables  order, 
p  =  IV7  remains  uniform  throughout  the  computational  domain. 


(b)  (d)  (f) 


Figure  3.1:  Recovery  of  fi(x)  (top)  and  the  approximation  error  (bottom)  from  their  N  =  128- 
mode  spectral  projections.  The  filter  (3.1)  was  of  orders  IV1/4  (a)-(b),  N 1//2  (c)-(d),  and 
max(2,  d( Nd(x ))4/2) 


(b)  (d)  (f) 


Figure  3.2:  Recovery  of  $2{x)  (top)  and  the  approximation  error  (bottom)  from  their  N  =  128- 
mode  spectral  projections.  The  filter  (3.1)  was  of  orders  IV1/4  (a)-(b),  IV1/2  (c)-(d),  and 
max(2,  ^(Nd(x))1/2)  (e)-(f). 
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(b) 


(d) 


(f) 


Figure  3.3:  Error  plots  for  the  recovery  of  fi(x)  (top)  and  f2(x)  (bottom)  from  their  N  = 
32,64,andl28  mode  spectral  projections.  The  filter  (3.1)  was  of  orders  A1/4  (a)-(b),  JV1//2  (c)- 
(d),  and  max(2 ,  ^(Nd(x))1^2)  (e)-(f). 


4.  Summary 

The  analysis  presented  here  quantitatively  resolves  the  classical  methodology  that  for  improved 
accuracy  low  order  filters  should  be  used  near  discontinuities  and  high  order  filters  away.  The 
optimal  adaptive  filters  presented  here  retain  the  traditional  robustness  associated  with  low  order 
filtering,  yet  achieve  a  significant  increase  in  accuracy  with  minimal  increase  to  computational  cost. 
Combined  with  the  automated  edge  detection  methods,  [GeTaOO],  adaptive  order  filtering  is  a  black 
box  procedure  for  the  exponentially  accurate  reconstruction  of  a  piecewise  smooth  function  from 
its  spectral  information. 
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